Spin gradient thermometry for ultracold atoms in optical lattices 
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We demonstrate spin gradient thermometry, a new general method of measuring the temperature 
of ultracold atoms in optical lattices. We realize a mixture of spins separated by a magnetic field 
gradient. Measurement of the width of the transition layer between the two spin domains serves 
as a new method of thermometry which is observed to work over a broad range of lattice depths 
and temperatures, including in the Mott insulator regime. We demonstrate the thermometry using 
ultracold rubidium atoms, and suggest that interesting spin physics can be reahzed in this system. 
The lowest measured temperature is 1 nK, indicating that the system has reached the quantum 
regime, where insulating shells are separated by superfluid layers. 

PACS numbers: 37.10. Jk, 03.75. Mn, 75.10.Jm, 05.30.Jp 



Ultracold atoms trapped in optical lattices represent a 
new frontier for the investigation of many-body physics 
[U [2] . The existence of novel physics at decreasing en- 
ergy scales drives the quest for lower temperatures in 
the atomic Mott insulator. Insulating Mott shells form 
at a temperature T ~ 0.2C/, where U is the interaction 
energy. At the lower temperature T ~ zJ, where J is 
the tunneling amplitude and z is the number of nearest- 
neighbors, the conducting layers become superfluid and 
the system enters a quantum insulator state |3j. At the 
even colder temperature scale T ^ J"^ /U , superexchange- 
stabilized phases can exist in the two-component Mott 
insulator; this is the regime of quantum magnetism 0]. 
Various proposals O |^ have focused on the realization 
of quantum spin Hamiltonians in this regime. Detection 
of superexchange-driven phase transitions in these sys- 
tems remains a major goal of ultracold atomic physics. 
Perhaps the most important barrier to experimental de- 
tection of such a phase transition is the requirement of 
temperatures well below 1 nK 4J. Additional cooling 
methods [THTU] will be needed to reach this very inter- 
esting temperature scale. However, it is clear that to 
assess current methods and to validate future cooling 
techniques, low-temperature thermometry of the Mott 
insulator is needed. 

Thermometry of systems in the Mott insulating state 
has remained a challenge [3l [TTkl4) . In this paper, we 
discuss and demonstrate a simple and direct method of 
thermometry using a magnetic field gradient which works 
in the two-component Mott insulator. 

The theory behind this method of thermometry is 
straightforward. The system under consideration is an 
ensemble of atoms in a mixture of two hyperfine states 
loaded into a three-dimensional optical lattice in the pres- 
ence of a weak magnetic field gradient. The two states 
have different magnetic moments, and are thus pulled 
towards opposite sides of the trapped sample by the gra- 
dient. At zero temperature, the spins will segregate com- 
pletely, and a sharp domain wall will exist between the 
two spin domains (a small width due to superexchange 



coupling is typically negligible). This system has the 
same bulk physics as the single-component Mott insu- 
lator, but includes additional degrees of freedom in the 
form of spin excitations in the domain wall. At finite 
temperature, spin excitations will increase the width of 
the domain wall. This width will depend in a simple way 
on the field gradient, the differential Zeenian shift, and 
the temperature, and can thus be used as a thermometer. 

For an incoherent mixture of two spins, the par- 
tition function for an individual lattice site can 
be approximately factorized as Z = Z^Zq^ where 
= I]<T6^P(~/5A*<T'B(a;)), /? is l/ksT, fi^ is the mag- 
netic moment of the spin cr, B(a;) is the spatially varying 
magnetic field, and Zq is the partition function of the 
particle- hole degrees of freedom (for which see [3j). This 
approximation is generally valid for the case of one atom 
per lattice site; for occupation number n > 1, it is valid 
when the mean of the intra-spin interaction energies Ua- 
is equal to the inter-spin interaction energy U-f^, which 
is a good approximation in *^Rb [15]. Since the total 
magnetization is fixed, the average value of the magnetic 
field is cancelled by the corresponding Lagrange multi- 
plier; we include this in the definition of B(a;). We are 
free to treat the two states as having pseudospin -|-1 and 
— 1; making that identification, the mean spin (s) as a 
function of position, gradient strength, and temperature 
has the simple form 

(s) =tanh(-/3- A/x-B(2;)/2), (1) 

where A/x is the difference between the magnetic mo- 
ments of the two states. A fit of the measured spin dis- 
tribution with a function of this form will give the tem- 
perature of the system. When the Zeeman shift due to 
the magnetic field gradient is a linear function of position, 
imaging of the spin distribution essentially corresponds 
to direct imaging of the Boltzmann distribution. 

The apparatus used to produce ultracold *^Rb atoms 
is described in Ref. il6 . After cooling, approximately 
10^ atoms are held in a far-red-detuned crossed optical 
dipole trap with trap frequencies between 100 and 200 
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FIG. 1: Images used for spin gradient thermometry. Data 
on the left were taken at a lower optical trap power than 
data on the right. Panels a and b are images of the spin 
distribution. Panels c and d show the mean spin versus x 
position. The fit to c gives a temperature of 52 nK; the fit 
to d gives a temperature of 296 nK. The inset of a shows the 
axes referred to in the text. The bar in b is a size scale. 



Hz. A three-dimensional cubic optical lattice, formed 
by three retroreflcctcd beams each of radius ~150/im, 
overlaps the trapping region. Since spin gradient ther- 
mometry does not depend on the number of atoms per 
lattice site n, we perform measurements at a range of n 
values between 1 and 4. The trapping and lattice beams 
are all derived from one fiber laser, with a wavelength A 
of 1064 nm. Magnetic field gradients up to a few G/cm 
can be applied with external coils, and calibrated using 
Stern-Gerlach separation of the different spin states after 
release from the trap. The gradient is applied along the x 
direction, which is the weakest axis of the crossed dipole 
trap. Absorptive imaging of the atoms is performed with 
a camera pointing down along the vertical z axis. 

The sequence of steps used to measure temperature 
is as follows. First, a sample of ^^Rb atoms in the 
\F = l,mF — — I) state is prepared by evaporation in the 
optical trap. Here F and mp are the quantum numbers 
for the total spin and its projection on the z axis, re- 
spectively. The atoms are then placed into a mixture of 
the |1, — 1) and |2, —2) states by a nonadiabatic magnetic 
field sweep through the microwave transition between the 
two states. This pair of states was chosen in order to 
avoid spin-exchange collisions. A magnetic field gradient 
of 2 G/cm is applied along the weak axis of the trap and 
results in additional evaporation, which is intended to re- 
move the entropy created by the state preparation |17j . 
At this point, the field gradient is changed to the value 
to be used for measurement; lower gradients are used 
for lower-temperature measurements to keep the domain 
wall width larger than the imaging resolution. The opti- 
cal lattice is then adiabatically ramped up, typically to 
a depth of 14.5 Er, where Er = h? /2mX^ is the recoil 
energy and m is the atomic mass. The transition to the 
Mott insulator occurs at 13.5 Er. At this point, the 
spin structure depends on the temperature as discussed 
above. 



There are several ways to measure the resulting spin 
distribution. One way is to first take an image of the 
F — 2 atoms in the 14.5 Er lattice, then in a second run 
to illuminate the atoms with an optical repumper beam 
resonant with the F = 1 to i^' = 2 transition for a few 
fis prior to imaging. This method gives an image of all 
atoms and an image of just the F = 2 atoms; appropriate 
subtraction can provide the spin distribution. It is possi- 
ble to determine the temperature from a single image of 
one spin, but the data in this paper were all taken using 
pairs of images to guard against systematic errors. 

The temperature can then be measured by fitting the 
spin distribution to the hyperbolic tangent form. The re- 
sulting thermometer has high dynamic range and variable 
sensitivity, works at all accessible temperatures of inter- 
est, and requires only the simplest fitting procedures. 
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FIG. 2: Independence of the measured temperature on the 
applied field gradient. The inverse of the width of the spin 
profile is plotted as a function of magnetic field gradient for 
two data sets at two different temperatures. For constant 
temperature, a linear curve is expected. The width is defined 
as the distance from the center to the position where the mean 
spin is 1/2. The solid (dashed) line assumes a temperature of 
123 nK (7 nK) and perfect imaging. The measured width of 
the colder data set saturates at high gradient because of finite 
imaging resolution. The dotted line assumes a temperature 
of 7 nK and an imaging resolution of 4 jj,m. 

Figure [T] shows data of the type used for spin gradi- 
ent thermometry. An image of the total atom density 
and an image of the spin density are obtained as dis- 
cussed above. Both images are then integrated along 
the y direction, which is transverse to the gradient. The 
spin distribution is then fit by a function of the form 
p{x) X tanh(|/?/iB^^^x), where p{x) is the total den- 
sity distribution. The only free parameters in this fit 
are a horizontal and vertical offset and the temperature 
T = l/kBp. 

Figures [2] and [3] show the results of this thermometry 
on ultracold ®^Rb atoms in an optical lattice. Figure [2] 
shows the linear scaling of the inverse width of the do- 
main wall as the magnetic field gradient is varied while 
holding the temperature constant. For widths larger than 
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the optical resolution, the scaling is as predicted by Eq. [T] 
The two data sets plotted in Fig. [2] were taken at two dif- 
ferent temperatures: 7 nK and 123 nK, according to the 
best-fit theoretical lines. Finite optical resolution or mo- 
tion of the atoms during imaging will blur the measured 
spin profile and result in an overestimate of the domain 
wall width at high gradients. This effect was modeled by 
applying a gaussian blur of radius 4 fim. to the theoreti- 
cal 7 nK spin profile at various gradients. The resulting 
curve, plotted as a dash-dot line in Fig. [2j reproduces 
the saturation of measured width observed in the exper- 
imental data. The effect of finite resolution is always to 
overestimate the temperature. 

Figure |3] shows the measured temperature plotted as 
a function of the power in the dipole trapping beam 
which confines the atoms in the direction of the mag- 
netic field gradient (the x direction). Higher powers in 
this beam lead to less effective evaporation, and thus 
higher final temperatures. As a check of the new method 
of thermometry, Fig. [3] also presents an analysis of the 
same data using an existing method of thermometry, 
based on measurement of the in-trap width of the atomic 
cloud along the direction perpendicular to the gradient. 
This second method is based on the well-known rela- 
tion (7^ = UbT /muP' ^ where a is the 1/e^ half-width 
of the atomic cloud and uj is the trap frequency in the 
direction along which the width is measured [H]. The 
width is determined by a fit to the wings of the trapped 
cloud. Trap width thermometry is based on a non- 
interacting approximation, and will fail at temperatures 
less than U when the system starts to become incom- 
pressible. As in Ref. [T2], all points on this plot are in 
the high-temperature single-band regime (T is less than 
the bandgap but greater than the bandwidth). For the 
temperatures plotted in Fig. [3j the agreement between 
the two methods is reasonably good, and gives confidence 
in the use of spin gradient thermometry in regions of pa- 
rameter space where no other thermometer exists. 

The large dynamic range of spin gradient thermome- 
try is evident in Fig. |3] Thermometry can be performed 
at temperatures so high that no condensate exists before 
lattice rampup. The lowest temperature we have mea- 
sured was achieved by using the new thermometry as a 
feedback signal, enabling adjustment of experimental pa- 
rameters for optimization of the final temperature in the 
Mott insulator. This method allowed us to achieve a mea- 
sured temperature as low as 1 nK. At the lattice depth 
used here, U is 37 nK, and zJ is 6 nK. The measured 
temperature is thus well below — zJ, the predicted 
critical temperature for the superfluid layer between the 
n — 1 and n = Mott domains. According to the treat- 
ment of Ref. , at 1 nK the system should be well inside 
the quantum regime, with concentric quantum insulator 
shells separated by superfluid layers. This represents the 
first direct demonstration that this temperature regime 
has been achieved in the Mott insulator. 
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FIG. 3: Validation of spin gradient thermometry. Comparison 
of two measured temperatures versus final power in one of 
the optical trapping beams. Squares represent the results of 
in-trap cloud width thermometry, and circles represent the 
results of spin gradient thermometry (see text for details). 
Error bars represent estimated uncertainties. The dashed line 
is a linear fit to the spin gradient thermometry data. The 
closeness of this fit suggests that the temperature reached is 
proportional to the trap depth. 



At a given value of the magnetic field gradient, very 
low temperatures will result in a width of the transition 
region smaller than the imaging optics can resolve (see 
Fig. [2]). However, the width can be increased by decreas- 
ing the magnetic field gradient. The lowest measurable 
temperature will then depend on the minimum achiev- 
able gradient as well as the optical resolution, which are 
technical rather than fundamental limitations. In our 
apparatus, background gradients with all coils turned off 
are of order 10~^ G/cm, which, given our imaging resolu- 
tion of a few /xm, would in principle allow measurement 
of temperatures down to ^ 50 pK or the superexchange 
scale, whichever is higher. 

It is instructive to compare the useful range of this new 
method of thermometry with that of existing methods. 
To facilitate meaningful comparison with non-lattice- 
based methods, we discuss the range of entropy per par- 
ticle S/Nks at which a given thermometer works, rather 
than the range of temperature. Condensate fraction ther- 
mometry works for 0.35 < S/Nks < 3.5, where the lower 
limit is set by the difficulty of detecting a thermal frac- 
tion less than 10%, and the upper limit is set by dis- 
appearance of the condensate. Thermometry based on 
the thermal cloud size has a similar lower bound, but 
extends to arbitrarily high values of S/Nkg. Quantita- 
tive thermometry based on the visibility of interference 
peaks upon release from the lattice requires state-of-the- 
art QMC calculations fitted to the data. This technique 
was recently used to measure temperatures as low as 
0.08J7 in the superfluid phase near the Mott insulator 
transition [18] . This method cannot be applied deep in 
the Mott insulating state PTI . Measurement of the width 
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of the conducting layers between the Mott shells is the 
only previously proposed method which works directly in 
the Mott insulating state [3l[4l[T9]. However, this method 
requires tomographic techniques, and the useful range of 
entropy is rather narrow: 0.4 < S/Nks < hi(2), where 
the upper limit is set by the melting of the Mott shells, 
and the lower limit is an estimate based on typical trap- 
ping parameters and optical resolution. Counting only 
spin excitations, the range of spin entropy per particle 
at which spin gradient thermometry works in our sys- 
tem is 0.1 < Sa/NkB < ln(2), where the lower limit is 
a function of optical resolution and sample size and the 
upper limit corresponds to the point at which the domain 
wall becomes as wide as the sample. It is important to 
note that spin gradient thermometry can work even if 
the entropy of the particle-hole excitations lies outside of 
this range in either direction. For example, spin gradi- 
ent thermometry can work at arbitrarily high values of 
the total entropy per particle S/Nks, assuming the field 
gradient is increased to the point where S^/Nks < ln(2). 

The method of thermometry presented here works be- 
cause the two-component Mott insulator in a field gra- 
dient has a spectrum of soft and easily measurable spin 
excitations. The wide dynamic range of this method is 
a result of the fact that, in contrast to the gapped spec- 
trum of the bulk one-component Mott insulator, the en- 
ergy of the spin excitations can be tuned by adjusting 
the strength of the magnetic field gradient. The addition 
of a field gradient and a second spin component docs not 
change the bulk properties of the Mott insulator and can 
be regarded as "attaching" a general thermometer to the 
first component. 

The two component Mott insulator in a field gradient 
is a rich system which can provide experimental access to 
novel spin physics as well as thermometry. In the work 
presented here, we have always allowed the spin distribu- 
tion to equilibrate in the gradient before ramping up the 
optical lattice. However, changing the gradient after the 
atoms were already loaded into the lattice should open 
up several interesting scientific opportunities, in which 
the gradient is used to manipulate or perturb the atoms 
rather than as a diagnostic tool. If, for example, the gra- 
dient were suddenly changed after lattice rampup, one 
could probe non-equilibrium spin dynamics in a many- 
body quantum system. If the gradient were instead low- 
ered adiabatically after rampup, adiabatic cooling of the 
Mott insulator could potentially be performed which, in 
contrast to |20| . would not involve spin- flip collisions. 

In conclusion, we have proposed and demonstrated a 



new method of thermometry for ultracold atoms in op- 
tical lattices. We have used the new method to measure 
temperatures in the Mott insulator as low as 1 nK. This 
temperature is to the best of our knowledge the lowest 
ever measured in a lattice, and it indicates that the sys- 
tem is deep in the quantum Mott regime. 
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